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I. INTRODUCTION 

We study theoretically the resistive state in short su- 
perconducting microbridges (with length L less than the 
coherence length ^(T)) at temperatures close to the crit- 
ical one. This subject was widely discussed and stud- 
ied in the 70-ies of last century (for review see and 
has regained renewed interest recently (see for example 
Refs. @|)- Interesting phenomena occurring in such a 
system, and which are still not completely understood, 
are: i) the hysteresis of current-voltage (IV) characteris- 
tics at relatively low temperatures [l| and ii) the 'foot'- 
like (sometimes also called 'shoulder'-like - see Ref. [1]) 
feature in the IV characteristics observed experimentally 
mainly in tin microbridges at low voltages and at tem- 
peratures close to Tc P-f?! (see also [ij). The hystere- 
sis is usually explained by Joule heating while for the 
'foot'-like structure several theories were proposed (8l-[l0| 
that are based on the idea that the quasiparticle distri- 
bution function is out of equilibrium (overcooled) as a 
consequence of the variation in time of the magnitude of 
the superconducting order parameter A = lAje**^ in the 
superconducting microbridge. The motivation for those 
theories comes from the fact that the energy of the quasi- 
particles depends on | A| and when the characteristic time 
scale for the variation of |A| is smaller than the inelastic 
relaxation time of the quasiparticle distribution func- 
tion /(e) the occupancy of the states with energy e may 
differ from the equilibrium one (for a detailed discussion 
see the book of Tinkham ) . 

To simplify the analytical treatment of the problem 
the authors of Refs. assumed that in the dynamic 
(resistive) state |A| varies as fast as S(j) and to find the 
coordinate and time dependence of \A\{x,t) they solved 
the stationary Ginzburg-Landau equation with a time 
dependent 6(j){t). Furthermore they assumed that: l)the 
relaxation term in the kinetic equations can be neglected 
when the period of oscillations of the order parameter is 



much smaller than Tin, 2) the spectral functions depend 
only on the local magnitude of the order parameter |A| 
(so called local approach). We argue that for a realistic 
inelastic relaxation time Tin, and even for a short micro- 
bridge L < ^(T) the nonequilibrium contributions to /(e) 
strongly affects the dynamics of | A| and results in a larger 
time scale for |A| than for Scf). Moreover the relaxation 
term in the kinetic equations plays a very important role 
at any voltage and therefore cannot be omitted. Taken 
together with the dependence of the spectral functions 
on both |A| and Q they provide an averaged 'heating' 
of the quasiparticles, instead of 'cooling' d, Q. Some 
kind of 'cooling' of quasiparticles at high voltages can be 
found only when one takes into account additional terms 
in the kinetic equations which couple the longitudinal /l 
(odd in energy) and transverse fx (even in energy) parts 
of 2/(e) = (1 — /^(e) — /T(e)) due to the finite spectral 
supercurrent (in previous works [1, these terms were 
omitted). This 'cooling' is not effective at low voltages 
and 'heating' together with the large time scale of the 
variation of | A| result in a hysteresis of the IV character- 
istics for relatively large Tin- We should stress that the 
time averaged 'heating' of quasiparticles is driven by os- 
cillations of |A| in the superconducting microbridge and 
not by Joule dissipation ^ j • E (j is the current density 
and E is the electric field). 

The paper is organized as follows. In section II we dis- 
cuss the theoretical model. In section III we present and 
discuss our results. In section IV we present our conclu- 
sions and discuss the possible origin of the experimentally 
found 'foot'-like structure in the IV characteristics. 



II. MODEL 

To simulate the resistive state in a short superconduct- 
ing microbridge we use the kinetic equations derived in 
Refs. [l^ - [l5l | for 'dirty' superconductors near Tc 
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here Q = {d<j}/dx — 2eA/c) is a quantity which is pro- 
portional to the superfluid velocity {vs = HQ/m), (p is 
an electrostatic potential, = II — fl ^^'^ /l(^) — 
tanh(e/2fcBT). iVi, N2, R2 are the spectral functions 
which should be found from the Usadel equation for the 
normal a(e) = cos0 = Ni{e) + iRi{e) and anomalous 
/3i = I3e"l>, P2 = Pe-"^ (/3(e) = sin 9 = N2{e) + zi?2(e)) 
Green functions 

hD^ + {{2ie — ^)-ni:>g2 cose) sin e+2|A| cose = 0, 

dX Tin 

(2) 

Eqs. l(a,b)are coupled through the finite spectral 
supercurrent [H [ill = i?e(/3iV^2 - ^2V^i)/2 = 
2A^2^2Q- Below we show that the coupling terms in 
Eqs. l(a,b) strongly infiuence S/l and the dynamics of 
the order parameter which is described by the modified 
time-dependent Ginzburg-Landau equation 
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The current and the electrostatic potential in the sam- 
ple can be found using the following equations 
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where (t„ is the normal state conductivity. In metals 
we have for the charge density p ~ and the condition 
divj = is satisfied due to Eqs. (lb, 3, 5). 

In the derivation of Eqs. l(a,b) it was assumed that 
deviations from equilibrium are small 5/l,/t ^ fl- 
it allowed one to linearize the collision integral due to 
electron-phonon collisions and to write it in the relax- 
ation time approximation. Furthermore, it was assumed 
that the inelastic relaxation time due to electron-electron 
interactions is much larger than due to electron-phonon 
hence one can neglect the corresponding collision in- 
tegral. 



where = nHD/SkBTc and A^^ = 

87r^(fcBTc)^/7C(3) are the zero temperature Ginzburg- 
Landau coherence length and the corresponding order 
parameter. Nonequilibrium parts of the quasiparticle 
distribution function enters Eq. (3) via the potentials 
$1 = - R2SfLde/\A\ and $2 = - N2fTde/\A\. 
When SJl is negative the potential $1 is positive and 
vice versus. In some respect, from the structure of Eq. 
(3) it follows, that one may introduce an 'effective' tem- 
perature for the quasiparticles Teff{x, t) = r-|-<I>i(a;, t)Tc 
and thus the positive/negative sign of $i(a;,i) means 
local 'heating'/'cooling' of quasiparticles. We should 
stress that we use the term 'effective' temperature only 
in order to give a simple physical interpretation of our 
numerical results and describe the integral effect of 
the nonequilibrium distribution /(e) which enters the 
equation for the order parameter via the potentials 
<^i{x,t) and $2(0;, i). Please note that the resulting 
nonequilibrium /(e) cannot be viewed as a Fermi-Dirac 
function with effective temperature Te//(x,t) and 
electrostatic potential ip{x,t). 
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FIG. 1: Schematic illustration of the model system with the 
geometric parameters. 

To model the superconducting microbridge connected 
to superconducting bulk leads in the current carrying 
regime we consider the system shown in Fig. 1. It 
consists of normal bulk reservoirs where the external 
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voltage Vext is applied and the constriction is modelled 
as a film with variable cross-section ^2 3> Si (here 
•Si (2) = Wi^2}di(2) is the cross-section, 1^1(2) is the width 
and di^2 is the thickness of the film (in this way we model 
variable thickness microbridge) . The narrowest part of 
the film is called the microbridge and the wide parts are 
assigned as being the leads. The wide part of the film of 
length 2L2 is in the normal state and the rest of the film 
is in the superconducting state (for relatively small Vext)- 
In our calculations we chose L2/L1 ^ 1 and the length 
Li was taken large enough to neglect nonequilibrium ef- 
fects from the NS boundaries located at a; = zt{Li + L/2) 
on the transport properties of the microbridge. The large 
normal part of the film is needed because in the bound- 
ary conditions for /l and /t enters the voltage V (see 
for example Ref. [l6|) and almost all the voltage drop 
occurs in the normal part of the film, effectively leading 
to a state with applied constant current even when the 
microbridge transits to the resistive state. Due to the 
large ratio S2/S1 » 1 the transport current does not de- 
stroy superconductivity in the wide part of the film and 
we mainly study the resistive state in the microbridge. 



In our calculations we assume that the lateral size of 
the film is much smaller than the London penetration 
depth A (or effective penetration depth = A^/c?i(2)) 
and therefore we may neglect screening effects. Further- 
more, we suppose that the current density is distributed 
uniformly over the width of the film. This assumption is 
strictly speaking not valid in the region where the cur- 
rent transits from the wide part to the narrow one (and 
vice versus). But because the current density is much 
smaller in the wide part than in the narrow part we may 
neglect such type of effects (at least for variable thickness 
microbridges) and consider variables which are averaged 
over the width and thickness of the film. Because the 
length of the microbridge is much larger than its width 
we may assume an uniform current distribution over the 
width of the microbridge. Thus our system is quasi ID 
and we only should take into account the continuity of 
A, 0, /t, SfL and conservation of the superconducting 
(Eq. 6(a) below), normal (Eq. 6(d)), energy (Eq. 6(c)) 
and spectral (Eq. 6(b)) currents at the points x = ±L/2 
through the boundary conditions 



S2dA/dx\±L/2 = SidA/dx\±L/2 (6a) 

S2de/dx\±L/2 = Side/dx\±L/2 (6b) 

S2D{{N! - Rl)d6fL/dx + jjT)\±L/2 - SiD{Nf - Rl)dSh/dx + jjT)\±L/2 (6c) 

S2D{{N^ + N^)dfT/dx + J,dfL)\±L/2 = SiD{{N^ + Rl)dfT/dx + J,6fL)\±L/2 (6d) 



It is easy to find the solution to Eqs. (1-5) in the 
normal region and therefore we only need to solve them in 
the superconducting region (in the interval —Li — L/2 < 
X < Li + L/2) with the following boundary conditions: 

^l±(Li+L/2) = 0, 6|±(Li+L/2) = 0, 5/l|±(Li+L/2) = 0, 

/t|±(Li+l/2) = ~eip± ■ dfl/dE, where ip± = TVext ± 
j±L2/(Jn and the current density j± in the points x — 
±{Li + L/2) could be found from Eq. (4). 

In our numerical calculations we use dimensionless 
units. The order parameter is scaled by Aq (Aq = 
l.lQkBTc ~ O.STAci is the zero temperature order pa- 
rameter value in the weak coupling limit), distance is 
in units of the zero temperature coherence length — 
^JKD/Kq ~ 1.2^(3i, time in units of = fi/ Aq and 
temperature in units of the critical temperature T^. The 
current is scaled in units of jo = Aocr„/(^oe), the super- 
fiuid velocity in units of Qo — ^c/2e^o a-nd the electro- 
static potential is in units of </7o = Ag/e. It is useful 
to introduce the dimensionless inelastic relaxation time 
^in = Tin/tQ which is the main control parameter in the 
model described by Eqs. l(a,b). 

We used the implicit method for the numerical solution 
of Eqs. (1-3). The coordinate step of the grid was equal 
to ^0 (which is much smaller than ^(T) for the consid- 



ered temperature interval 0.92 < T/T^ < 0.99) and the 
time step varied from 0.5 to up to 2 depending on the 
temperature. In our numerical procedure we use an even 
number of grid points and consequently the center of the 
microbridge where the order parameter goes to zero (the 
phase slip center) is situated between grid points. There- 
fore any quantity in the phase slip center (PSC) is in fact 
calculated at a distance 6x = ^o/2 ^ £,{T) from the PSC. 

For the geometrical parameters of the film we used 
the following values: 52/51 = 10, L2 = 500fo and Li 
was varied from 15 ^0 (at T = 0.92Tc) up to 45 ^0 (at 



T = 0.99Tc). In our calculations we used fi, 



1000 



which covers typical values for many low temperature 
superconductors (for example 15| in Nb fi„ 10^ and 
in Al fin — 10'^ ). In the wide part of the film we used 
fin = 5 which allows us to decrease Lin and the length 
Li in order to optimize the calculation time. 

To find the current-voltage characteristics of the su- 
perconducting microbridge we applied a large voltage 
iVext to the normal reservoirs which induces a large 
current (/ > Ic) resulting in a resistive state in the mi- 
crobridge. Than we decrease Vext in a step-wise man- 
ner and we find the time averaged difference of the elec- 
trostatic potentials between the ends of the microbridge 
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FIG. 2: (color online) Time dependence of the current, ampli- 
tude of the order parameter in the center of the microbridge 
and voltage drop over the microbridge at fixed external volt- 
age Vext- The results were found from numerical solution of 
Eqs. (1-5) in the presence of coupling terms in Eqs. l(a,b). 
The geometrical parameters of the system are: W2/W1 — 10, 
L2 = 500^0, Li = 22^0 and L = 7^0- 

{V = Tp{—L/2) — Tp{L/2)) as function of Vext- In a sim- 
ilar way it is easy to find the critical voltage V^^t (or 
critical current Ic) of the superconducting microbridge 
when the superconducting state becomes unstable (by 
increasing Vext from small values). As a result we find 
the dependence V {Vext /Vext) which practically coincides 
with the dependence V{I/Ic)- In Fig. 2 we plot the 
time-dependence of the current in the microbridge which 
illustrates that indeed / varies weakly in time and the 
constant Vext induces an almost constant / in our model 
system even when the microbridge is in the resistive state. 



III. RESULTS 

In Fig. 3 we show the current-voltage characteristics 
of a short microbridge with length L — 7^o at T = 0.96Tc 
calculated for different values of the inelastic relaxation 
time Tin in the regime of decreasing applied current (the 
length and temperature were chosen close to the param- 
eters of the experiment of Ref. Q and for tin fm ^ 200). 
Firstly, we should note the strong infiuence of the cou- 
pling terms in Eqs. l(a,b) on the IV characteristics: volt- 
age increases (decreases) with increasing rj„ at I > Ic in 
case of the absence (presence) of coupling terms. Sec- 
ondly, in both cases the IV curves are hyperbolic-like 
with no 'foot '-like feature and with hysteresis for rela- 
tively large Tin (the transition to the superconducting 
state occurs at the retrapping current Ir which could be 
smaller than the critical current Ic for the transition from 
the superconducting to the resistive state). 

We plot in Fig. 4 the energy dependence of (5/l aver- 
aged over one oscillation period T|a| taken in the phase 




FIG. 3: (color online). Current-voltage characteristics (in the 
regime of decreasing current) for a short microbridge with 
length L — 7^o at T = 0.96Tc calculated for different values 
of the inelastic relaxation time in the absence (a) and the 
presence (b) of coupling terms in Eqs. l(a,b). Solid curve is 
the IV characteristic of a short microbridge in the absence of 
nonequilibrium effects (see for example [J]). 



slip center. The effect of the coupling terms are clearly 
visible by comparing Fig. 4(a) and 4(b). The (5/^ is 
on the average negative in case the coupling terms are 
absent and positive when they are included. 

Let us now discuss the origin of the sign of (5/^ for 
the different energies. First we consider the case when 
the coupling terms are neglected. Let us simplify Eq. 
1(a) to 

dt nn Ni de dt ' ^ ' 

where we omit the diffusive and coupling terms. If we 
average Eq. (7) over r|A| then we obtain 

If i?2 and Ni are solely a function of |A| then 5/^ — 0. 
But in general the spectral functions R2 and V^i are a 
function of two variables |A|, \Q\ (at fixed x and e - see 
Eq. (2)). To get insight how the ratio R2/N1 changes 
with varying jQj at fixed |A| one may solve Eq. (2) with 
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e > \A\{x,t){e < |A|(x,t)) than in the decreasing re- 
gion. It results in Sfj^ < in the phase shp center for 

e>|A| and <5/i > for e < |A| (see Fig. 4(a)). 
The time averaged potential $i is positive in the micro- 
bridge (see inset in Fig. 4(b)) which implies 'heating' 
of quasiparticles and suppression of the superconducting 
properties. The larger Tin the stronger the deviation form 
equilibrium (see Eq. (8) and Figs. 4(a)) and the larger 
the hysteresis (see Fig. 1(a)). 
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FIG. 4: (color online). Typical energy dependence of the time 
averaged S/l in the center of the short microbridge (i.e. in 
the phase slip center) in the absence (a) and the presence (b) 
of the coupling terms in Eqs. l(a,b). In the inset of Fig. 4(b) 
we show the distribution of the time averaged "I>i and |A 
along the microbridge and the leads (schematically shown). 



zero second derivative and find that the ratio R2/N1 de- 
creases when e > I A I and increases when e < | A | for large 

\Q\- 

In Fig. 5 we present the time dependence of |A|, Q and 
$1 in the absence of the coupling terms in Eqs. l(a,b). 
First of all we should note that the amplitude of the os- 
cillations of |A| is smaller than Aiead- The reason for 
this effect is the large characteristic time scale for the 
variation of |A| in comparison with the one of 6(f>. In- 
deed it is known that nonequilibrium effects may con- 
siderably slow down the dynamics of |A| [l7l - [l9| . When 
A I increases (decreases) its value is smaller (larger) than 
one could expect from the static dependence \A\stat{S(t>) 
(see for example Eq. (1.3) in Q) for the given value of 
6(j)(t) ^ Q(t)L. For example when 6(f> reaches zero (which 
corresponds to the moment in time when Q = after the 
phase slip event in Fig. 5) the order parameter is still 
much smaller than Aiead- 

Due to the time delay in the variation of |A| the su- 
pervelocity Q ^ 6(1)/ L is different for the same values 
of |A| taken at different times during the period T|^| 
(see dashed lines in Fig. 5). Therefore the ratio R2/N1 
in the increasing region of |A| is larger (smaller) for 



FIG. 5: (color online). Time dependence of the order pa- 
rameter |A|, the supervelocity Q and the potential $1 in the 
center of the short microbridge (parameters are the same as 
in Fig. 4) at / = 1.2/c and Tin = 250 as found from a nu- 
merical solution of Eqs. (1-3) without coupling terms. The 
horizontal dashed line marks |A| = Aiead/2 to illustrate the 
corresponding different Q values at two moments in time. 

Consider now the effect of the coupling terms in Eqs. 
l(a,b). The term D\/{jjT) = 2D\/ [N2R2Q fr) m 
Eq. 1(a) may be considered as an additional source 
of nonequilibrium which is nonzero in the energy in- 
terval 5e ~ near local |A|(a;) (in the spatially uni- 
form case with Q — Q and Tin = 00 one has N2R2 ~ 
|A|(j;)(5(e — |A|(a;)) - see for example Ref. [§]). From 
Eq. 1(b) we may roughly estimate /t ^ —etpdf^/de and 
the sign of the coupling term DS/^jefx) is defined mainly 
by the product Q • E (for discussion of the effect of the 
coupling terms see also Ref. [19])- 

Our calculations show (see Fig. 6) that near the super- 
conducting leads the product Q ■ E is mainly (in time) 
positive and in the phase slip center it can be both pos- 
itive and negative (due to the sign change of Q - see 
Fig. 5). It results in positive 6fL at e ~ Ai^ad - see 
Fig. 4(b) (at these energies the ends of the microbridge 
are determinative because near the leads |A| has a low 
oscillation amplitude and Q • E > 0). At lower energies 
e ^ \A\pg^ < Aiead the main source of nonequilibrium 
is situated in the center of the microbridge where both 
Q-E and d\A\/dt changes sign during the oscillations 
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FIG. 6: (color online). Time dependence of the product Q ■ E 
in the phase slip center and near the ends of the short super- 
conducting microbridge during T^a\ as found from a numerical 
solution of Eqs. (1-3) including the coupling terms. 



PSC 

and S/l is negative except at low energies e ^ |A| 

(see Fig. 4(b)). The net effect is a negative time averaged 
$1 (see inset in Fig. 4(b)) and therefore an enhance- 
ment of the superconducting properties which explains 
the decrease of the voltage at fixed current I ^ Ic with 
increase of Ti„ (see Fig. 3(b)). But at low voltages the 
effect of the coupling terms (^ QE) becomes smaller and 
the role of 'heating' increases. This subsequently leads 
to an increase of the voltage at fixed current / ~ /c with 
increasing Ti„ and to an hysteresis of the IV curves (but 
which is smaller than in the case without coupling terms 
- compare Figs. 5(a) and 5(b)). 



IV. DISCUSSION 

The origin of the here predicted time averaged 'heat- 
ing' of quasiparticles is different from Joule dissipation. 
Indeed, the source of nonequilibrium ~ d\A\/dt changes 
sign during oscillations of the order parameter while the 
term I ■ V is always positive. Time averaged 'heating' 
appears only due to the dependence of the spectral func- 
tions both on |A| and Q and due to the difference in 
time variation of |A| and 5(f>. Energy for 'heating' of the 
quasiparticles (and 'cooling' due to the term S/{jefT) in 
the kinetic equation) comes from the energy / • V deliv- 
ered by the external source of the current (this energy 
goes also for example to the 'heating' of phonons due to 
presence of the relaxation term in the kinetic equations). 

Our result is opposite to the one of Refs. [1, Q where a 
time averaged 'cooling' of quasiparticles and an enhance- 
ment of superconductivity was found at high voltages 
(when T\A\ ^ '''in) in the absence of the coupling terms 
in Eqs. l(a,b). The differences with our calculations 



are that we take into account: i) the backaction of the 
nonequilibrium quasiparticle distribution function on the 
dynamics of |A|; and ii) the dependence of the spectral 
functions A^i, N2 and i?2 on both |AJ and \Q\. There is 
also one more difference with Refs. 0,0 - we did not ne- 
glect the relaxation term in Eq. (7) at high voltages (see 
Eq. (3) in Ref. @ and Eq. (16) in Ref. [i)^By omitting 
this term one may obtain 'cooling' with Sfj^ > if one 
chooses the initial condition S/l = when |A| = Aiead 
in the center of the microbridge (and when the spectral 
functions i?2and A^i depend only on | A | as was assumed 
in Refs. @, Q). But this choice of the initial condition 
is not obvious. By choosing a different initial condition 
(large negative SJl at |A| = Ajead) one may obtain 'heat- 
ing' in this case too. Only the presence of the relaxation 
term in Eq. (7) resolves this problem, because it de- 
fines undoubtedly the time averaged value of S/l (see 
Eq. (8)). Note that due to model assumptions for spec- 
tral functions of Refs. [1, ^ 5f has to be equal to zero 
at any voltage if one does not neglect the relaxation term. 
From a physical point of view the relaxation term in Eq. 
(7) tells us that the accumulation of the 'heat'/'cold' (af- 
ter sudden change of the current for example) occurs on 
a time scale which is proportional to Tin ^ T\a\ and one 
should take it into account to find the true value of 5f 
and T|^| corresponding to the given value of the current. 
It provides a strong dependence of the IV characteristics 
and 5f on even at high voltages (see Figs. 3(a) and 
4(a)) while in Refs. no dependence on Tin was found 
in this limit. 

According to our numerical calculations the IV char- 
acteristics of short 'dirty' superconducting microbridges 
are hyperbolic-like and hysteretic (at relatively large Tin 
and/or low temperatures) in the temperature interval 
0.92 < T/Tc < 1 (calculated and not presented here 
the IV characteristics at r/rc=0.92-0.99 are qualitatively 
similar to those of Figs. 3(a,b)). We believe that the 
experimentally found 'foot'-like features in the IV char- 
acteristics of short superconducting microbridges is con- 
nected with the geometry of the studied microbridges - 
they were short and wide (in experiments [sKtI the width 
W ~ 2L). For such a geometry the current density dis- 
tribution over the width of the microbridge (constriction) 
should be nonuniform with sharp peaks near the edges 
of the microbridge In Ref. [2l[ it was shown that a 
nonuniform current density distribution in the supercon- 
ductor may lead to slow vortex motion at low currents 
and phase slip lines at high currents [2l| and to an IV 
characteristic (see for example Fig. 4 in Ref. 21|) which 
resembles the experimental results of Refs. p-H]- More- 
over the sharp transition from low to high voltages oc- 
curs at Vc which is inversely proportional to Tin |2l| - the 
same result was found in the experiment on tin super- 
conducting microbridges Q- Therefore, we predict that 
the IV characteristics of short {L < £,{T)) superconduct- 
ing microbridges should change from hyperbolic-like (in 
case when W <^ L) to the one with a well-pronounced 
'foot'-like structure at low voltages when W > L. Good 
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candidates are Sn, In, Pb and Al with relatively large 
coherence length (^(T) > 500nm at T ~ TJ. We ex- 
pect that our results will be tested in the near future 
because recently it has become possible to fabricate mi- 
crobridges with length L < 500 nm and width W > 50 
nm (see for example Ref. 0,0]) which satisfy the condi- 
tion W <s:L< ^(T) at 0.9 < T/T^ < 1. 
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